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Astrophysical observations have put unprecedentedly tight constraints on cosmological theories. 

The ACDM model, mathematically simple and fits observational data-sets well, is preferred for 
explaining the behavior of universe. But many basic features of the dark sectors are still unknown, 
which leaves rooms for various nonstandard cosmological hypotheses. As the pressure of cosmological 
constant dark energy is unvarying, ignoring contributions from radiation and curvature terms at low 
redshift, the effective pressure keeps constant. In this paper, we propose two parametric models 
for non-constant effective pressure in order to study the tiny deviation from ACDM at low redshift. 

We recover our phenomenological models in the scenarios of quintessence and phantom fields, and 
explore the behavior of scalar field and potential. We constrain our model parameters with SNe la 
and BAO observations, and detect subtle hints of Ude < — 1 from the data fitting results of both 
models, which indicates possibly a phantom dark energy scenario at present. 


I. INTRODUCTION 


Since the discovery of current acceleration of our uni¬ 
verse expansion in 1998, maybe the greatest mystery in 
cosmology is the deceptive nature of the dark energy. Re¬ 
cent observational results [l| have put tight constraints 
on the properties of dark energy, but there is still no 
theoretical or observational indication to pin down its 
nature. On one hand, although the simple cosmological 
constant A can accommodate the accelerating expansion, 
it encounters two serious problems. The first one is the 
fine tuning problem: the measured energy of the vacuum 
so much smaller than the estimated value 
which is the famous 120-orders-of-magnitude discrepancy 
that makes the vacuum explanation suspecious. The sec¬ 
ond one is why the dominance of the cosmological con¬ 
stant over the matter component at the present epoch. 
These two basic problems prompt us to propose some al¬ 
ternatives, which include an evolving scalar field called 
quintessen ce [2 ~B, noncanonical scalar field (such as K- 
essence [9l-[llI7 phantom 0, i modified grav¬ 

ity 0j S i 191423j| , coupled dark energy 0, [23, or de¬ 
caying dark energy [26j | and so on. On the other hand, 
as we know the equation of state (EoS) parameter of the 
cosmological constant is precisely aide = — 1- Recent ob¬ 
servations show that the EoS parameter of modeled dark 
energy is Wde = —1.006 ± 0.045, which slightly favours 
ojde < —1- Anyhow, the small deviations from the cos¬ 
mological constant A allow one to consider models with 
Wde 7^—1- So one can make efforts to construct new mod¬ 
els to explain the deviations which may be detectable at 
the precision of current and future observations. 
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Parameterization is an useful tool towards a more com¬ 
plete characterization of dark energy modelling and has 
been routinely employed to analyze datasets. Most pa- 
rameterizations for dark energy models involve the EoS 
parameter ujde for the dark energy behavior. Several well- 
known parameterizations for the EoS of dark energy have 
been proposed so far. We can write parameterizations in 
polynomial form uJde{z) = u!nXn{z) generally, where 

n—0 

the expansions can be given by the following, (i)Redshift; 
Xn{z) = z", (ii)Scale factor: Xniz) = (1-^)” = (ify)", 
(iii)Logarithmic: Xn(z) = [ln(l-|-z)]". Parameterization 
(i) was proposed by Huterer and Turner}^ and Weller 
and Albrecht with n < 1. Parameterization (ii) with 
n < 1 was introduced by Chevalier, Polarski and Lin¬ 
der the famous Chevallier-Polarski-Linder(CPL) 

parameterization, ujde = wo + d;i(l —a) = ujq+ujij^ be¬ 
haves as uJde -T ujQ+uJi for z —>■ oo and ojde —t wq for z —>■ 
0. A more general form with ojde = wq + was 

later proposed by Jassal, Bagla and Padmanabhanfs^. 
Parameterization (iii) with n < 1 was introduced by 
Efstathiou[3lj- In recent years, some new parameteri¬ 
zations have been proposed, such as using Pade parame¬ 
terizations for the EoS of dark energy}^, namely uJde = 
T+tlfi-a? >and i.de = It is worth mentioning 

that Sen proposed a parameterization for the pressure of 

dark energy model [sl, 113, = —-Po + Pi(l — a) H-, 

in order to study small deviations from the cosmologi¬ 
cal constant. Different from parameterizations which fo¬ 
cused on the EoS of dark energy mentioned above, in this 
paper we aim to make parameterizations for the relation 
between redshift and effective pressure of all energy com¬ 
ponents in the universe. In the following we proposed two 
parametric models for the effective pressure in order to 
explore late-stage evolution of the universe. 

This paper is organized as follows: in section II, we 
propose two new parametric models for the effective pres- 
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sure: P{z) = Pa + PbZ and P{z) = Pc + In 

Section III, we relate our parametric models with the 
quintessence and phantom scalar fields, and the behavior 
of field and potential is then explored. In Section IV, we 
constrain our model parameters with SNe la and BAO 
observations. In Section V, We end with discussions and 
conclusions. 


II. TWO PARAMETRIC MODELS 

The Friedmann equations, equation of energy conser¬ 
vation and equation of state constitute a close system to 
describe the background evolution of the universe. A sub¬ 
stitute from EoS to a relation between effective pressure 
P and redshift z is also feasible, as equation P = P{z) is 
not linearly dependent on the Friedmann equation and 
equation of energy conservation. Also, the EoS can be 
recovered by inserting P — z relation into equation of 
energy conservation 

p + 3H{P + p) = 0, (1) 

and integrating out the expression of p. For example, 
the effective pressure for ACDM at late stage is nearly 
constant, say Pq; accordingly, we can obtain from Eq. o 
that 


p{a) = -Po + Ca-^ (2) 

where C is an integration constant, and the two terms at 
the right side represent contributions from cosmological 
constant and matter respectively. 

This is just an example for P parameterization; gen¬ 
erally, we can have more complicated P — z relations. 
As P — z relation is equivalent to EoS, a parameteriza¬ 
tion on the effective pressure is equivalent to that on the 
EoS parameter ujde- Since ujde is the exponential of some 
component in EoS, uJde paramerterization requests a pre¬ 
supposition of the components in EoS; i.e., the physical 
mechanism of the possible deviation from ACDM has to 
be dictated although We make parameterizations merely 
because we actually do not know the concrete mechanism 
behind the accelerative expansion. To illustrate, a devia¬ 
tion of ACDM might come from the evolution of the EOS 
of the cosmological constant term, while an additional 
component can also result in same deviation. However, 
a parameterization on the effective pressure just circum¬ 
vents this issue, and do not require any knowledge of the 
concrete physical mechanism. We are able to directly 
study the deviation from constant P — z relation without 
prejudice to a presupposition. 


A. Model 1 

In this subsection, we propose a model which reads: 


(3) 


where Pa and Pb are free parameters. 

For scale factor a and redshift z, we have 


a = 


do 

1 + z 


1 

1 + z' 


(4) 


where oq = 1 corresponds to the value today. Substitute 
Eqs. ©and (jH) to Eq. (P), the total energy density can 
be integrated as 

P(a) = -(Pa - Pb) - ^Pb^~^ + ( 5 ) 


where Ci is an integration constant. If we set poto be the 
energy density today, the integration constant is then 
Cl = Po + Pa + ^Pb- In Eq. (P), we can interpret the 
inversely cubic term Cia~^ as dust matter and the con¬ 
stant term —{Pa — Pb) as the cosmological constant in 
ACDM. Term does not appear in the ACDM 

model, whose physical nature will be explored in next 
section. 

For convenience in date htting, we introduce some di¬ 
mensionless parameters. First, we define dimensionless 
density and pressure as 


P* 

P* 


P__^ 

7o~M' 

p 

Po 


( 6 ) 

(7) 


The expressions of the total density Eq. and total 
pressure Eq. can be rewritten as: 


p*{a) = -{P: - p;) - + C*a-^ (8) 

P*{a) = {p:-P;) + p:a-\ (9) 

where P; ^ P; ^ ^ and Q ^ ^ = I + P* + \P;. 

Redefined two new parameters a = — (D* — P^) and 
(3 = —^Pb, then: 


p*{a) = a + 13a ^ + {l — a — j3)a (10) 

P*{a) = (II) 

As we know, the dimensionless Hubble parameter is 


E{z) ^ (12) 

tlo 

Compare Eq. (EH) with Eq. ®, we obtain, 

E{a) = p*{a)^. (13) 


Then, for model I, we define 


Hi 

O 2 


a 

f3a~^ 

UmOd 

E2 


(14) 

(15) 

(16) 


where flmo = 1 — a — P, hence Di -|- H 2 -I- flm = 1- 


P{z) = Pa + PbZ, 
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B. Model 2 

We propose another parameterization as 

P(z) = Pc + ^, (17) 

i -h Z 

where Pc and Pd are free parameters. Inserting Eqs. o 
and (El) into Eq. du, we obtain the total energy density 
for model 2 , 


p{a) =-Pc-^Pda +C 2 a (18) 

where C 2 is an integration constant. Set the present en¬ 
ergy density as po, then C 2 = Po + Pc + \Pd- Still, we can 
find term C 2 a~^ corresponding to dust matter, and term 
—PcCorresponding to the cosmological constant. The dif¬ 
ference between model 2 and model 1 rests on the rest 
term, — |Pha“^ for model 2 whereas —jPda for model 
1 .Their physical nature will be studied in next section. 

Like model 1, we need to introduce new model parame¬ 
ters in model 2. With Eqs. dH]) and ([7]), we can obtain the 
expressions of total density and total pressure for model 
2 : 


p*{a) = -P*-^^P*a + C*^a-\ 

(19) 

p*{a) = p; + p>. 

( 20 ) 

where P* P*, ^ f and ^ ^ = 1 + P* 


Redefine two new parameters 7 = — P* and <5 = 
then: 

3 p* 

4^d ’ 

p*{a) = 7 - 1 -(5a- 1-(1 — 7 — ^)a“^. 

( 21 ) 

P*{a) = - 7 -^( 5 a. 

( 22 ) 

Also, we define for model 2 , 


II 

(23) 

0 - 

2 £,2’ 

(24) 

^ ^jYlQd 

(25) 


where flmo = 1 — 7 — <5, we have fli -f f 22 + Hm = 1- 


III. RELATION WITH SCALAR FIELDS 

Deviations from the ACDM in our models can be real¬ 
ized through different physical scenarios. Scalar fields are 
mainstream approaches to explain the acceleration of the 
universe expansion. In the scenarios of scalar fields, dark 
energy evolves with time. The scalar field dynamics has 
been studied by literature in great details (see Ref. UM) 
and there are lots of issues involved such as (i)choosing 


initial conditions for scalar field; (ii) choosing potential 
with solid theoretical motivation; (iii) the existence of 
the tracker field and so on. Generally, evolution of scalar 
field is studied over the cosmic history, and once the pa¬ 
rameters of scalar field models are set they determine the 
entire cosmological evolution. So a more detailed anal¬ 
ysis would involve studying scalar field dynamics over 
cosmic history, and then comparing its evolution with 
that of pressure parametrization model at low redshift. 
In this paper, we will merely compare pressure and en¬ 
ergy density of field with that of a model of pressure 
parametrization at low redshift and study the behavior 
of field and potential. Physical realization of parameteri- 
zations through scalar fields means adjusting behavior of 
scalar fields to dark energy term occurred in parametric 
model. Specihcally speaking, we make two equations 


Peff — Pscalar field: (^(^) 

Peff Pm — Pscalar field: (^7) 

as mathematical definition of realization. 

In this section, we will take “quintessence” and “phan¬ 
tom” as two examples to realize our models. 

Quintessence: “Quintessence” denotes a canonical 
scalar field (f) with a potential Vi {(p) that does nor interact 
with all the other components except standard gravity, 
whose EoS parameter ujde > ~1- The quintessence is 
described by action 

5 = j d^X^[^R + ^^]+SM: ( 28 ) 

■ 5^0 = (</.), ( 29 ) 

where = 87 rG, R is the Ricci scalar and Sm is the 
action of matter. The variation of the action Eq. (1^^ 
with respect to (j) gives 

E/(?i) =0, (30) 


where Vi(()') is the potential of the quintessence field, 
the prime denotes the derivative with respect to (j). In a 
FLRW background, the energy density pde and the pres¬ 
sure Pde of the quintessence field are 


Then the EoS 


Pde — 


(31) 

II 


(32) 

^de — 

- Vi{4>) 

(33) 


Phantom: Minimally coupled phantom model is also 
a possible realization, whose EoS parameter uide < ~1- 
The action of the phantom field minimally coupled to 
gravity and matter sources is 


5 = j d^X^[^R + ^^]+SM: (34) 

= ]^9^''d^(j)d:,(t)-V2{(l)): 


( 35 ) 
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whose variation with respect to </> gives 
0 + = 0 , 


we choose condition 4>a=i = Mp, where Mp is the re¬ 
duced Planck mass. The Friedmann equation can then 
(36) be rewritten as 


where V 2 { 4 >) is the potential of the phantom field, and 
the prime denotes the derivative with respect to (j). The 
energy density and pressure of the phantom are given 
by(assuming flat FRW metric) 

Pde = — 2 *^^ 

Pde = 

The EoS of the phantom field is then 

-lJ^-V2i^) 

^de = -T-:-■ 

-^<P^+V2{cl>) 

So UJde < -1 for < V 2 {(l)). 


A. Model 1 

The EoS of the scalar fields for model 1 reads 

Pscalarfield i/?(l + z) 

“ -— —i H- r~a7T~, —T- 

Pscalarfield ^ “t” /t(l T 

Note that in above equation, there will be a singularity 
when z = — 1. In this paper we only consider the 

universe at low redshift, so we need not to worry about 
that situation. Besides in section IV data fitting will 
support our argument. 

In the quintessence scenario, assuming the cosmic com¬ 
ponents consist of matter and quintessence, comparing 
Eq. ^ cLuci Eq. (|32p with Eq. (|3]) S/iid. Eq. @, we have 

- {Pa - Pb) - ^Pba-^ = i<^' + Ei(0), (41) 

Pa-Pb + Pba-^ = i<^^-Fi(0). (42) 

Simplify the above two equations, compare to 
Eqs. © 0 (0) dSl) Cni) (HH), replace model parame¬ 
ters {Pa , Pb) with redefined parameters {a , /3), then we 
obtain 


(37) 

(38) 

(39) 


= ^Po/3a \ (43) 

5 

Vi((/>) = poa + -pQl3a~^. (44) 

From Eq. (|4^ . it is easy to find that /3 > 0 in the sce¬ 
nario of quintessence. By Eq. (H51) and Eq. (Hill , one can 

construct the kinetic energy and potential Vi{cj)) of 
the quintessence field with parameters {a , /3) of model 1. 
In order to solve the above two equations, following [3^ . 




1 


(45) 


Considering the dark energy domination at present 
epoch in the universe, with the density parameter in dark 
energy ~ 0.7, we define Vq = po = ‘iMpH^. Simplify 
Eq. (H51) and Eq. (HU) , we have 


d(j) 

da 


—— = zLMpt 




a 


+ jda ^ + {I — a — P)a ^ 


^(46) 


Vi{^) = 


Ho(a + -/3a ). 


(47) 


The symbol “±” in Eq. (l46l) corresponds to two solu¬ 
tions. Consider a = 0.7,/3 = 0.05 for numerically solving 
the above two equations, the solutions are represented 
in Fig. o and Fig. 0 respectively. From Fig. o , we 
can find </> increases with a, the potential decreases with 
the increasing </>, Eq. (03), implying that the potential 
will reach the minimum value Fl(()>) = Voa in the fu¬ 
ture. From Fig. 0, we can see that (j) decreases with 
a, and the potential decreases with the decreasing </>, 
Eq. (1471) implies that the potential will reach the min¬ 
imum Vi{(j)) = Voa in the future. By Eq. (ITTl) and 
Eq. (flSl) . we can obtain the expression of density pa¬ 
rameter for model I: 


= III + ' ■ (48) 

The evolution of density parameter in the scenario 
of quintessence is plotted in Fig. 0. From Fig. © , we 
can see that until low redshift the energy density in the 
quintessence field becomes cosmologically dominant. Fi¬ 
nally, the field comes to rest at the minimum of the po¬ 
tential Vi{(j)) = Voa, and the universe eventually settles 
into a de Sitter phase (see Eq. (gni)). 

In the case of phantom scenario, assuming the cosmo¬ 
logical components consist of matter and phantom, com¬ 
paring Eq. (1371) and Eq. (I38p with Eq. 0 and Eq. ©, 
then we have 


- {Pa - Pb) - ^na-1 = + F2(<^), (49) 

Pa-Pb + Pba-^ = V2(0). (50) 

Replace model parameters {Pa , Pb) with redefined pa¬ 
rameters (a , /3), we have 


= -^Po/3a \ 

(51) 

5 

V 2 ((/') = poa + -pol3a~^. 

(52) 
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FIG. 1: The solution of Eq. (H51l and Eq. (H71) corresponding 
to a plus sign in Eq. (1461) . Field (p as a, function of a depicted 
in the top panel, potential Vi as a function of (j> depicted in 
the bottom panel. The arrow indicates evolutional direction 
of potential with respect to time. We consider values a = 
0.7,/3 = 0.05. 


FIG. 2: The solution of Eq. (1461) and Eq. (1471) corresponding 
to a minus sign in Eq. (M . Field (p as a function of a depicted 
in the top panel, potential Vi as a function of (p depicted in 
the bottom panel. The arrow indicates evolutional direction 
of potential with respect to time. We consider values a = 
0.7,13 = 0.05. 


From Eq. dsn, it is easy to find that in the scenario of 
phantom, /3 < 0. By Eq. (EU and Eq. (l5^ . one can 
construct the kinetic energy and potential V 2 {(p) of 
the phantom field with model parameter (a , /?). Eq. (ISTl) 
and Eq. dSH) can be rewritten as 


d(p 

da 

V2{4>) 


dzA/p4 


; + /3a ^ + (I — a — j3)c 


-3 


a-^(53) 


Vb(a + -/3a ^). 
6 


(54) 


Consider a = 0.7,/3 = —0.05 for numerically solving 
the above two equations, the two solutions are repre¬ 
sented in Fig. (HI) and Fig. ([S]) respectively. From Fig. (jH), 
we can find ^ increases with a, and the potential increases 
with the increasing (p, Eq. (iMl) implies that the poten¬ 
tial will reach the maximum value V 2 {(p) = Vqcx in the 
future. From Fig. (p decreases with a, and the po¬ 
tential increases with the decreasing cp, in the future the 
potential will reach the maximum value V 2 {(p) = Voa. In 
Fig. dSl), we plot the evolution of density parameter in 
the scenario of phantom. Notice that the energy density 



Z 


FIG. 3: Evolution of the density parameters in the 
quintessence field and matter (flm) for model 1. 0,,/, 

is indicated by solid line, and flm is indicated by dashed line. 
We consider values a = 0.7, /3 = 0.05. 


in the phantom field becomes cosmologically dominant 
only in the recent past. In the future, the field comes to 
rest at the maximum of the potential and the universe 
eventually settles into a de Sitter phase. 
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FIG. 4: The solution of Eq. ((53ll and Eq. (l54)) corresponding 
to a plus sign in Eq. (I53II . Field as a function of a depicted 
in the top panel, potential Vi as a function of (j> depicted in 
the bottom panel. The arrow indicates evolutional direction 
of potential with respect to time. We consider values a = 
0.7,/3 = -0.05. 


FIG. 5: The solution of Eq. (1^ and Eq. (ImI) corresponding 
to a minus sign in Eq. (O. Field (j) as a. function of a depicted 
in the top panel, potential Vi as a function of (j> depicted in 
the bottom panel. The arrow indicates evolutional direction 
of potential with respect to time. We consider values a = 
0.7,/3 = -0.05. 


B. Model 2 

Write down the EoS of the scalar fields for model 2: 


^scalar field 

Wde = -^- = -l- .An ^-l ^ 

Pscalarfield T “t" o(l Z) 

It is obvious that only when parameters ( 7 , 5) have op¬ 
posite signs, there will be a singularity coming out when 
z = —1 — In section IV, data htting results will show 
that such a singularity would not appear at low shift. 

In the quintessence scenario, compare Eq. m and 
Eq. (15^ with Eq. (IT^ and Eq. (IT71) . we can obtain 

(56) 

Pc + Pdo = El(<(>). (57) 

Simplify the above two equations, referring to 
Eqs. ® © (HU) (EHl) (ED) (EH) replace model parame¬ 
ters (Pc , Pd) with redefined parameters (7 , S), then: 



Z 


FIG. 6 : Evolution of the density parameters in the phantom 
field (B<^) and matter (fim) for model 1 . is indicated by 
solid line, and flm is indicated by dashed line. We consider 
values a = 0.7, (3 = —0.05. 


= -^Po<5a, 

(58) 

7 

= Po7+g/0o^a- 

(59) 
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From Eq. (l58l) . it is easy to find that in the scenario 
of quintessence <5 < 0. By Eq. (155)) and Eq. (1551) . the ki¬ 
netic energy and potential Vi((/)) of the quintessence 
field are constructed with parameters (7 , S) of model 2 . 
Simplify these two equations, we have 


d(j) 

da 


zLMpi 


-S 


7 -f (5a -f (1 — 7 — S)a 


-3 


7 (60) 


Vii(l>) = Vo{-f+'^Sa), ( 61 ) 

where Vq = po = iMpH^. Choose parameters 7 = 
0.7, (5 = —0.05 for numerically solving the above two 
equations, the two solutions are represented in Fig. dZD 
and Fig. dS]) respectively. From Fig. ( 0 , we can find (j) 
increases with a, and the potential decreases with the in¬ 
creasing (j). From Fig. (0, we can find (j) decreases with a, 
and the potential decreases with the decreasing (j). Notice 
that since (5 < 0 in the scenario of quintessence, according 
to Eqs. ( 0 ) (EJ) (113, the Friedmann equation is written 
as = 3]^Po[7 + ^a -I- (1 — 7 — 5)a“^], which will not 
hold when the scale factor a is very large. Nevertheless 
at low redshift the relation is still feasible. 

By Eq. (l23l) and Eq. (1241) , we can obtain the expression 
of density parameter for model 2: 


= + = (62) 

the evolution curve of density parameter in the sce¬ 
nario of quintessence has been plotted in Eig. ([9]), from 
which we see the quintessence field begins to dominate 
at low redshift. 

In order to realize model 2 in a phantom scenario, com¬ 
paring Eq. (1571) and Eq. (1551) with Eq. (1T51) and Eq. (1T71) , 
we can obtain 



a 



<t> 


FIG. 7: The solution of Eq. (I 6 OI 1 and Eq. (l6T1) corresponding 
to a plus sign in Eq. (1601) . Field (j) as a function of a depicted 
in the top panel, potential Vi as a function of (j> depicted in 
the bottom panel. The arrow indicates evolutional direction 
of potential with respect to time. We consider values 7 = 
0.7, (5 = -0.05. 


d(j) 

da 


zLMpt 


7 -I- (5a -I- (I — 7 — (5)c 


-3 


(67) 


-P,-^P,a = -i(^2+E2((^), (63) 

Pc + Pda = (64) 

Simplify and replace model parameters {Pc , Pd) with 
redefined parameters (7 , i5), we have 

= ^Po<5a, (65) 

V2(</>) = Pol+'^Po^a. (66) 

6 

From Eq. (1551) and Eq. (1551) . it is easy to find that 
in the scenario of phantom (5 > 0. By the above two 
equations, one can construct the kinetic energy and 
potential V 2 ((/>) of the phantom field with parameters (7 
, i5) of model 2. Eq. (1551) and Eq. (1551) can be rewritten 
as 


= Vb (7 + ^5a). ( 68 ) 

6 

Choose parameters 7 = 0.7, S = 0.05 for numerically 
solving the above two equations, the two solutions are 
represented in Eig. m and Eig. dni) respectively. Erom 
Fig. (Uni), we can find (p increases with a, and the po¬ 
tential increases with the increasing cp. In Fig. dni), cp 
decreases with a, and the potential increases with the 
decreasing p. Notice that since (5 > 0 in the scenario 
of phantom, the Friedmann equation can be written as 
= ■^^p[)['y + 5a + {\ —7 — (5)a“^], H —>■ 00 as a —>■ 00 , 
which means there will be a “rip” in the future. 

In Fig. m, we plot the evolution curve of density 
parameter in the scenario of phantom. Note that the 
phantom becomes cosmologically dominant only in the 
recent past, finally the EoS parameter uJde is less than 
— I {uJde = ~|, see Eq. (1551) 1 and the universe eventually 
settles into a “rip”. 
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a a 




<l> 




FIG. 8 : The solution of Eq. (l60ll and Eq. (|6T1) corresponding 
to a minus sign in Eq. dSOl). Field as a function of a depicted 
in the top panel, potential Vi as a function of (j> depicted in 
the bottom panel. The arrow indicates evolutional direction 
of potential with respect to time. We consider values 7 = 
0.7, J = -0.05. 


FIG. 10: The solution of Eq. (IBTIi and Eq. corresponding 
to a plus sign in Eq. iszi). Field as a function of a depicted 
in the top panel, potential V 2 as a function of (j> depicted in 
the bottom panel. The arrow indicates evolutional direction 
of potential with respect to time. We consider values 7 = 
0.7,5 = 0.05. 



shift range z = [0.015,1.4]. To perform the chi-square 
statistics, the theoretical distance modulus is defined as 


fith{zi) = 5\ogiQDL{zi) + ^lo, (69) 


where ^0 = 42.39 —5 log^g h with h the Hubble parameter 
in units of lOOkm/s/Mpc, 


Dl = 



dz' 

E{z'-,ey 


(70) 


z 


FIG. 9: Evolution of the density parameters in the 
quintessence field (14,^) and matter (flm) for model 2 . 14,^ 
is indicated by solid line, and f4m is indicated by dashed line. 
We consider values 7 = 0.7, <5 = —0.05. 


is the Hubble-free luminosity distance in a spatially flat 
FRW universe, E{z-, 9) is the dimensionless Hubble pa¬ 
rameter, and 9 is model parameters. 

The corresponding Xsn function is calculated from 

2 [fj-obsizi) — fj,th{zi)y 

Xsn — / . 7 ^ I 

i=l t 


IV. ASTROPHYSICAL DATA CONSTRAINTS 

Type la Supernovae In this paper we use the 
Union2.1 SNe la data-sets without systematic errors for 
data fitting, which compiles 580 SNe la covering the red- 


where piobs{zi) and <Ji are the observed value and the 
corresponding Icr error of distance modulus for each su¬ 
pernova. The minimization with respect to po can be 
made trivially by expanding Xsn 

Xsn = -\- /igC, 


(72) 
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a 



<P 


580 


m =i: 


^obs{zi) - = 0 ) 


580 




(74) 

(75) 


Thus /j,o is minimized as /xq = ^ by calculating the 
following transformed ■ 

xl^[9)=A{e)-Z^. (76) 

Baryon Acoustic Oscillations The baryon acoustic 
oscillation (BAO) data-sets are listed in Table (P). We use 
the parameter A to measure the BAO peak in the distri¬ 
bution of SDSS luminous red galaxies. In the following 
A is defined as 

. /— - , ,_ir 1 dz' ,2 

A=V^oEizt) 3[-y^ ^]a, 

where zt = 0.35. The for BAO data is 

2 \Aobs{zi) — Ath{zi\9)Y 

XbAO - 2^ ^2 ■ 

i=l A 


(77) 


(78) 


FIG. 11: The solution of Eq. (l67ll and Eq. (IM)) corresponding 
to a minus sign in Eq. (EZl). Field </) as a function of a depicted 
in the top panel, potential Vi as a function of (j> depicted in 
the bottom panel. The arrow indicates evolutional direction 
of potential with respect to time. We consider values 7 = 
0.7, J = 0.05. 



Z 


FIG. 12: Evolution of the density parameters in the phantom 
field (fl</)) and matter (Hm) for model 2 . is indicated by 
solid line, and is indicated by dashed line. We consider 
values 7 = 0.7, S — 0.05. 


The total is given by 

= As AT + Xbao- (’i’9) 

The fitting results and corresponding reduced for 
model 1 and model 2 are listed in Table (P. The like¬ 
lihoods of parameter (a , /3) and (7 , S) are shown in 
Fig. (fO|) and Fig. (ITTl) . respectively. Besides, evolution of 
the EoS parameter ujde with respect to redshift z with la 
error propagation from data-fitting (Tab. P are shown 
in Fig. (IT5|) and Fig. (HU), respectively. 

TABLE I: 6 measurement points of the Baryon Acoustic Os¬ 
cillation Data-sets. 


redshift 

A aA 

Sample 

0.106 

0.526 0.028 

6 dFGS [37] 

0.20 

0.488 0.016 

SDSS [32 

0.35 

0.484 0.016 

SDSS [37] 

0.44 

0.474 0.034 WiggleZ 

0.6 

0.452 0.018 WiggleZ [32 

0.73 

0.424 0.021 WiggleZ [32 


V. CONCLUSION 


where 


580 


A(e) = S 


iMobs(zi) - fithizi; 0; fio = 0 )]^ 


(73) 


Since the observational confirmation on late-stage ac¬ 
celerative expansion of the universe many years ago, dif¬ 
ferent models have been proposed to explain its source, 
among which parameterization is a widely used scheme to 
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TABLE II: parameters of model 1 and model 2 estimated by 
SNe la and BAO data-sets with la errors, 
model 1 model 2 

xLn/d-o./. 564.045/(583) xLn/d.o.f. 564.098/(583) 
a 0.771 ±0.084 7 0.635 ±0.119 

P -0.058 ±0.084 5 0.079 ±0.120 



FIG. 13: Icr and 2a confidence ranges for parameter pair 
(ce, P) of model 1, constrained by SNe la and BAO data¬ 
sets. The dotted straight line(/3 = 0) corresponds to a ACDM 
model. The blue dotted line and red dotted line correspond 
to flmo = 0.26 and flmo = 0.32, respectively. 



r 


FIG. 14: la and 2a confidence ranges for parameter pair ( 7 , 5) 
of model 2, constrained by SNe la and BAO data-sets. The 
dotted straight line((5 = 0) corresponds to a ACDM model. 
The blue dotted line and red dotted line correspond to flmo = 
0.26 and flmo = 0.32, respectively. 


better characterize the dark energy with observational re¬ 
sults. In this paper, we studied two models parameteriz¬ 
ing the effective pressure at low redshift, P{z) = Pa + PbZ 
and P{z) = Pc + 

Deviations from the ACDM can be realized through 
different physical scenarios. Roughly speaking, there 
are two ways. One is to introduce some small but 
nonzero components besides the cosmological constant 



z 


FIG. 15: Evolution of the EoS parameter ujde as a function 
of the redshift z with la error propagation, constrained by 
SNe la and BAO data-sets for model 1. The solid line, the 
straight dotted line, and light blue region represent the best- 
fit, ojde = —l(ACDM), and la region, respectively. 



FIG. 16: Evolution of the EoS parameter iJde as a function 
of the redshift z with la error propagation, constrained by 
SNe la and BAO data-sets for model 2. The solid line, the 
straight dotted line, light blue region represent the best-ht, 
ujde = — l(ACDM), and la region, respectively. 


A, such as imperfect fluid cosmology [3§ - l^ and cosmic 
strings ; whereas the other is to assume the cos¬ 

mological constant A exactly zero and the dark energy 
characterized by scalar fields evolving with time. In this 
paper, we pick the second way. We presented two parain- 
eterizations in the scenarios of quintessence and phantom 
fields, and accordingly expressed the kinetic energy term 
and potential term V{4>) with model parameters (a , 
P) and (7 , (5) respectively. Then we reconstruct the den¬ 
sity parameter for quintessence and phantom evolv¬ 
ing with redshift. In order to obtain a better physical 
understanding of the field and potential, we numerically 
solved the field as a function of the scale factor a and the 
potential as a function of field (f). 

We constrained model parameters {a , /3) and (7 , 
i5) with SNe la and BAO data-sets. We reconstructed 
evolution of the EoS parameter oj^e in term with the 
redshift z. For model 1, the value for EoS parame¬ 
ter UdeO is -1.027'^:[j;‘Q^43 at present day; for model 2, 
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Wdeo = ~1-0371 q;q5q. These results show that model 1 
and model 2 both slightly indicate that the EoS param¬ 
eter of dark energy ojde < —1, which corresponds to a 
phantom dark energy scenario at present. Still, we can 
not rule out a quintessence dark energy scenario or a A 
dark energy scenario. 

Different parameterizations possess their own advan¬ 
tages in addressing some particular problems, but their 
validity may not be ensured when applied for explaining 
global evolution. For example, our two parameterizations 
on effective pressure can estimate the deviation from the 
prediction of standard model at low redshift with gen¬ 


erality that does not depend on the concrete physical 
mechanism behind. 
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